Optical linewidth of a low density Fermi-Dirac gas 
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We study propagation of light in a Fermi-Dirac gas at 
zero temperature. We analytically obtain the leading den- 
sity correction to the optical linewidth. This correction is a 
direct consequence of the quantum statistical correlations of 
atomic positions that modify the optical interactions between 
the atoms at small interatomic separations. The gas exhibits 
a dramatic line narrowing already at very low densities. 
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The observation of Bose-Einstein condensation in di- 
lute atomic vapors m has spurred much interest in ultra- 
cold atomic gases. Another evident milestone of atomic 
physics would be the cooling of a Fermi-Dirac (FD) gas 
to the quantum degenerate regime. So far all probing 
of atomic Bose-Einstein condensates has been done opti- 
cally, and obviously optical detection could also play an 
important role in the experimental studies of FD gases. 
Appropriately, theoretical studies of the FD gases are ex- 
periencing a renaissance &-@] . 

In this paper we investigate propagation of low- 
intensity light in a FD gas in the limit of low atom den- 
sity. We derive the leading quantum statistical correction 
to the standard column-density refractive index analyt- 
ically, by legitimately ignoring collective linewidths and 
line shifts generated in the processes in which a photon is 
repeatedly scattered between the same atoms. A fermion 
gas exhibits a striking line narrowing characteristic of the 
FD statistics, which behaves as p^^^ at low atom densi- 
ties. 

In the dipole approximation it is advantageous to 
transform the Hamiltonian into the length gauge by the 
Power-Zienau-WooUey transformation |0]. Then the pos- 
itive frequency component of the electric field E+ may be 
expressed P,B[ in terms of the positive frequency compo- 
nents of the driving electric displacement, DJ^, and of the 
source field radiated by atomic polarization, P"*", as 
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Here k — fJ/c, Vt is the frequency of the driving field, 
and the scalar constant k = D^ /heo is defined in terms 
of the reduced dipole moment matrix element V. The 
monochromatic dipole radiation kernel G(r) coincides 
with the corresponding classical expression ^^. In sec- 
ond quantization the polarization reads 
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Here ipg and V'e are the ground state and the excited state 
atom field operators in the Heisenberg picture, and dge 
is the dipole matrix element for the transition g — > e. For 
simplicity, we consider here two-level atoms with just a 
single ground state \g) and one excited state |e), using a 
constant real vector d (such that T) — |d|) as the dipole 
matrix element. 

While Eqs. @) describe the scattered light in a 
medium, in general, with a small atom-light detuning and 
for a dense atomic sample, there is no easy way to find 
the polarization P"*" (r) . By making a field theory version 
of the Born and Markov approximations, we have derived 
a hierarchy of equations of motion for correlation func- 
tions that contain one excited-atom field and one, three, 
five, etc., ground state atom fields, for the limit of low 
light intensity |g] . In the present case of two-level atoms 
the hierarchy reads 

P/(ri,---,r/_i;r;) = (iJ - 7)P;(ri, • • • , r;_i;r;) 
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where 7 = T)^ k^ / (GTrheo) denotes the spontaneous 
linewidth, and S is the atom-light detuning. We have 
defined a projection operator 
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whose purpose is to eliminate all but the two atomic 
states from consideration, and the correlation functions 

P;(ri,- ••,r;_i;r/) 

^ (V'J(ri) . . . ^t(r,_i)P+(r07/^,(r,_i) ■ • • ^^(ri)), (5a) 

p,(ri, . . . , r,) ^ (V;t(ri) • ■ • ^J(r0^g(rO • • ■ ^'^(ri)). (5b) 

The quantity P; reflects correlations between the dipole 
moment of one atom and the positions of ^ — 1 other 
atoms, and pi is simply the density correlation function 
for / ground state atoms. 

The terms in the sum on the right-hand side of Eqs. (H) 
represent processes in which the I atoms at ri, . . . , r; re- 
peatedly exchange photons. Such processes are the mi- 
croscopic mechanism for collective linewidths and line 



shifts. The integral stands for a process in which yet 
another atom shines its hght on the atom at r; . 

Due to the resulting divergent dipole-dipole interac- 
tions, all correlation functions P/ vanish whenever two 
position arguments are the same M . The Lorentz-Lorenz 
local-field correction follows mathematically from this 
observation. Moreover, without changing the outcome 
of the hierarchy, we may, and will, remove all contact 
interactions between different atoms in Eqs. (0) by intro- 
ducing the field propagator G' defined by 



G^(r) = Gy(r)+m<5,,<5(r)/3. 
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This definition indicates that the integral of G' over an 
infinitesimal volume enclosing the origin vanishes. 

The coupled theory for light and matter fields [Eqs. (|l|) 
and (|3|)] may be solved, in principle exactly, by means of 
stochastic simulations Q. This is because the correla- 
tion hierarchy (^) is the same as the hierarchy describ- 
ing the classical electrodynamics of charged harmonic os- 
cillators with the position correlations pi. By synthe- 
sizing a stochastic ensemble of samples of dipoles that 
have the position correlation functions pi and calculat- 
ing the ensemble-averaged response to classical light, we 
have a solution to Eqs. (||). Unfortunately, such simula- 
tions are demanding on computer time. The computa- 
tions of Ref. H were therefore performed within a one- 
dimensional (ID) model electrodynamics. While the pre- 
dictive power of ID electrodynamics may be questioned, 
the simulation results for a FD gas at T = show clear 
signatures of the quantum statistics: Even in the limit of 
zero density, the optical linewidth of the FD gas is only 
half of the resonance linewidth of an isolated atom M . 

The ID simulations have also allowed us to test pre- 
dictions of the density expansion introduced by Morice 
et al. [|ll| in their studies of the optical response of a 
quantum degenerate Bose-Einstein gas. At least in one 
dimension this expansion is in semi-quantitative agree- 
ment with numerical simulations, and in the low-density 
limit the agreement is excellent [H . With this in mind, we 
venture to use the approximation of Morice et al. |ll| ] to 
truncate the correlation hierarchy (Bl) also in the present 
three-dimensional case. 

We consider the steady-state solution of (H). The 
atoms are assumed to fill the half-infinite space z > 
with a constant density p. The incoming free field is 
written Di^(r) = Dp eexjp{ikz),a,nd we assume that 
e||d. The hierarchy of equations (y) is truncated by writ- 
ing(l| 

P3(ri,r2;r3)c^^^^^^i^P2(r2;r3). (7) 

P 

The pair correlation function is written in the form 

P2(ri,r2)-p2[l + </j(ri-r2)]. (8) 

By introducing the dimensionless quantities 
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where a denotes the dimensionless atomic polarizability, 
and the ansatz Pi(r) = Peexp(iA;'z) with Im(fc') > 0, 
we obtain the susceptibility of the sample as 
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Here we use the dimensionless integration variable f = 
kr. The quantity C, and hence also fc', have been forced 
to be independent of position by essentially ignoring the 
effects of the surface of the atomic sample [^ . 

The second term in the denominator in Eq. (HO) gives 
the Lorentz-Lorenz shift. In the absence of the C term 
the electric susceptibility is the standard column den- 
sity result augmented with a local-field correction. Atom 
statistics and atom-field collective effects are encapsu- 
lated in the integral C. The expansion of Ref. ||ll|] is 
such that the parameter C takes into account quantum 
statistical position correlations between any pair of atoms 
and the exchange of photons between them to arbitrary 
order, but ignores all repeated photon exchange involving 
more than two atoms. Correspondingly, one may expand 
the functions inside the integrals in Eq. (O) as power se- 
ries in aG', and interpret the nth order as a contribution 
in which a photon is radiated between a pair of atoms n 
times. 

In this paper we only consider the leading order in the 
modifications of the optical response of an atom due to 
the presence of the other atoms, and write 
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The expression (h2f) arises from processes in which any 
"probe" dipole is subject to the external driving field, 
and in addition to the primary radiation from the other 
dipoles. Collective linewidths and line shifts, processes 
that involve the repeated scattering of a photon between 
the same atoms, are ignored. For uncorrelated locations 
of the dipoles with t/j = 0, the effects of the primary ra- 
diation from the other dipoles on a probe dipole average 
to zero. However, for the FD statistics Eq. (O) gives a 
nontrivial result. This reflects the short-range ordering, 
within the correlation length, of the atoms in the gas. 

Fermions at T = fill the Fermi sphere fik = 8(/cf — 
|k|), with the Fermi wave number kp = {Q'k'^pY''^. In 
this case we may evaluate the pair correlation function 



[Eq. ^] in closed form. In the thermodynamic hmit the 
result is 
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We also set the polarization of the incoming light field to 
be in the x direction. The propagator in Eq. ( p^ ) then 
reads 
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After inserting Eqs. (|l^ and (^ into Eq. (|l|) we ob- 
tain a (complicated) analytical expression, whose density 
expansion reads 
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In our ID electrodynamics, the entire expression (nV) 
may be integrated analytically for a FD gas at T = [^ . 
It is then easy to see that the lowest-order density con- 
tribution is correctly introduced by the expansion (12). 
This is also the case in three dimensions, although the 
demonstration is more indirect. First, the terms in 
Eq. (|ll| ) that do not depend on ip{r) are linearly pro- 
portional to p. Second, expanding the contributions to 
Eq. ( [ll|) that do depend on ip into a series of aG', for 
orders beyond the one included in dl2) we find radial 
integrals of the form 

df r^ ^;—^ip{f) oc 0(p), 
dff^^ip{f)o^O{pY'^ . 



Here ^ are integers independent of p, and n > 3. Most 
of these integrals formally diverge at the origin, but in a 
manner that must eventually cancel to give a finite re- 
sult ([ll|). Besides, the divergences do not depend on den- 
sity. Omitting the divergences, the integrals scale with 
density as indicated. All told, Eq. (|lj) not only repre- 
sents the lowest-order correction to the optical properties 
in terms of the number of microscopic optical interaction 
processes between the atoms, but it also correctly gives 
the leading density correction ex p^/'^ to the optical re- 
sponse. 

We have plotted the linewidth and the line shift, in- 
cluding both the effect of the FD statistics and the 
Lorentz-Lorenz shift, 
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as a function of density in Fig. |] using the full form of 
C from of Eq. (|lj) . A zero-temperature FD gas exhibits 
a striking linewidth narrowing already at low densities. 
For p = 1.5 X 10"'^ the optical linewidth of the gas is 
r ~ 0.947, and with p = 1.5 x 10"^ we have T ~ 0.79 7. 
For the 767 nm transition in ''"K [g| the corresponding 
densities would be p ~ 8.2 x 10^^ cm^^ and p ~ 8.2 x 10^^ 
cm~'^. At p = 0.1 (p ~ 5.5 X 10^^ cm~^) the optical 
linewidth of a FD gas would be approximately half of 
the linewidth of an isolated atom. However, at the lat- 
ter density we may already have to consider collective 
linewidths and line shifts to obtain a reliable quantita- 
tive prediction, a task we do not undertake in the present 
paper. From Fig. |](b), the line shift is negative at low 
densities and completely vanishes at p ~ 0.03. At higher 
densities the line shift turns positive. 

The dramatic line narrowing may be attributed to the 
regular spacing between the atoms characteristic of the 
FD statistics. As discussed in Ref. Q, the mechanism is 
particularly transparent in one dimension. An alterna- 
tive description of the line narrowing may be obtained in 
the momentum representation. At T = the fermions 
fill the Fermi sphere. Due to the Pauli exclusion prin- 
ciple only either strict forward scattering or scattering 
events that take the recoiling atom out of the Fermi sea 
are allowed. The change of the wave vector of an atom 
upon scattering satisfies |Ak;| = 2ksm{9/2), where 6 is 
the scattering angle for photons. All atoms are scattered 
out of the Fermi sea if \Ak\ > 2fci?. Thus, we see that 
for the photon scattering angles 6 satisfying 



sin(6'/2) > kp/k = (67r2p)i/3 
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scattering is not inhibited by the FD statistics. On 
the other hand, for sin(0/2) < kp/k some recoil events 
would lead to an already occupied state in the Fermi sea, 
and are forbidden. The suppression of light scattering 
is strongest in the near-forward direction corresponding 
to small values of 9. When the density is increased, at 
p > l/(67r^) we have kp > k, and scattering is at least 
partially suppressed in all nonforward directions. Corre- 
spondingly, Fig. n shows that p = l/(67r^) ~ 0.017 is a 
relevant scale for the density. 

It is instructive to note the difference between differ- 
ent atom statistics. For the Bose-Einstein condensate 
the standard factorization of the correlation functions, 
pi — p'', corresponds to an uncorrelated atomic sample, 
and gives 1^ = 0. The leading correction to the standard 
column density linewidth therefore results from the co- 
operative optical effects, the collective optical linewidths 
and line shifts, and it is proportional to atom density, 
just as the Lorentz-Lorenz shift. FD statistics is differ- 
ent because the correlations have a length scale kp^ that 
enters the argument, and the length scale itself depends 
on density: kp^ oc p^^/'^. The result is that, at low den- 
sities, the effects of the FD statistics dominate over the 
Lorentz-Lorenz shift. A Maxwell-Boltzmann ideal gas 
has another nontrivial correlation function ip, but the 



length scale is determined by temperature and does not 
depend on density. The leading density correction to the 
optical response is then once more proportional to atom 
density. 
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FIG. 1. The optical (a) linewidth and (b) the line shift of a 
Fermi-Dirac gas as a function of the atomic density per cubic 
optical wave number of the driving light p = p/k^ . 

We have assumed a homogeneous FD gas in our anal- 
ysis. Now, a FD gas may be considered locally homoge- 
neous lH] when the length scale over which the density 
varies is much larger than the spatial correlation length. 
Given the length scale of a harmonic trap I = {fi/mujY/'^ 
and the correlation length 1/kp — (Gtt^p)"^/'^ from 
Eq. ([l3|), the criterion reads l^p ^^ 1. A simple dimen- 
sional argument shows that this is the same as requiring 
that the number of trapped atoms be much larger than 
one. Furthermore, if the size scale of the atomic sample 
is much larger than the wavelength of light, / ^ A, it 
is reasonable to expect that the refractive index as ap- 
propriate for the local density applies in the bulk of the 
gas. 

In conclusion, we have studied propagation of light in 
a FD gas. We have discussed quantum statistical correc- 
tions to the refractive index, and calculated the leading 
density correction to the standard column density sus- 
ceptibility. Already at low densities, fermions exhibit a 
dramatic narrowing of the resonance line. This might 
serve as a signature of quantum degeneracy in a cold FD 
gas. 
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